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ABSTRACT 

We acquired spectra of 141 H ii regions in ten late-type low surface brightness galaxies 
(LSBGs). The analysis of the chemical abundances obtained from the nebular emission 
lines shows that metallicity gradients are a common feature of LSBGs, contrary to 
previous claims concerning the absence of such gradients in this class of galaxies. 
The average slope, when expressed in units of the isophotal radius, is found to be 
significantly shallower in comparison to galaxies of high surface brightness. This result 
can be attributed to the reduced surface brightness range measured across their discs, 
when combined with a universal surface mass density-metallicity relation. With a 
similar argument we explain the common abundance gradient observed in high surface 
brightness galaxy (HSBG) discs and its approximate dispersion. This conclusion is 
reinforced by our result that LSBGs share the same common abundance gradient with 
HSBGs, when the slope is expressed in terms of the exponential disc scale length. 
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1 INTRODUCTION 

During the past few decades the investigation of the ra¬ 
dial chemical abundance gradients of spiral galaxies has 
provided critical constraints for our understanding of their 
evolutionary histories. The overall picture that emerges 
from a large body of observational work on the oxy¬ 
gen abundances (O/H) of extragalactic Hll regions, the 
standard tracers of metals in star-forming galaxies (see 
Pilyugin, Grebel & Kniazev 2014a for a recent compila¬ 
tion), is consistent with an inside-out growth of galac¬ 
tic discs, in which the abundance gradients reflect radial 
variations of star formation rates and gas infall timescales 
(Prantzos & Boissier 2000; Fu et al. 2013). 

Recent observational work on the chemical composi¬ 
tion of spiral galaxies has addressed a number of topics, 
in particular the use of abundance diagnostics that are 
complementary to Hll regions, such as blue supergiant 
stars (Kudritzki et al. 2008, 2012; Bresolin et al. 2009b) 
and planetary nebulae (Bresolin et al. 2010; Stasihska et al. 
2013), azimuthal variations (Li, Bresolin & Kennicutt 
2013; Sanchez et al. 2015), the flattening of the gradients 
at large galactocentric distances (Bresolin et al. 2009a; 
Bresolin, Kennicutt & Ryan-Weber 2012; Marino et al. 
2012) and in interacting systems (Kewley et al. 2010; 
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Rich et al. 2012), and the spatial distributions of metals in 
high-redshift galaxies (Yuan et al. 2011; Queyrel et al. 2012; 
Jones et al. 2013). Moreover, integral field spectroscopy is 
currently providing abundance gradients for large samples of 
nearby galaxies (Rosales-Ortega et al. 2010; Sanchez et al. 
2014), yielding high statistical significance to correlations 
with global galactic properties and scaling relations. 

Despite this progress, the abundance gradient prop¬ 
erties of low surface brightness galaxies (LSBGs), a 
significant component of the overall galaxy population (e.g. 
Impey & Bothun 1997), remain virtually unknown, largely 
as a result of observational challenges. Chemical abun¬ 
dances of late-type LSBGs have insofar been obtained either 
from emission line spectroscopy of the central portions 
of compact systems (Burkholder, Impey & Sprayberry 
2001; Bergmann, Jprgensen & Hill 2003; Liang et al. 
2010) or from individual Hll regions, but with 
insufficient numbers to derive reliable radial 
trends (McGaugh 1994; Rbnnback & Bergvall 1995; 
Kuzio de Naray, McGaugh & de Blok 2004). To date the 
only attempt to constrain the oxygen abundance gradient 
in LSBGs has been reported by de Blok & van der Hulst 
(1998), who, examining nebular abundances in three 
LSBGs, concluded that their radial distributions are fiat. 

The typical oxygen abundances of LSBGs span 
a wide range, from 12-|-log(0/H) < 8.4 (McGaugh 
1994; Ronnback & Bergvall 1995) up to super-solar val- 
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ues^ (Bergmann et al. 2003), but on average peak¬ 
ing below the characteristic abundances measured in 
high surface brightness galaxies (Liang et al. 2010). 
For example, Burkholder et al. (2001) found a mean 
value 12-f log(0/H) = 8.22 for 17 LSBGs, compared to 
12-I-log(0/H) = 8.47 for 49 high surface brightness galax¬ 
ies. 

The relatively low metallicities, together with small 
star formation rates and blue optical and near-infrared 
colors, corroborate the notion that LSBGs are slowly- 
evolving systems (de Blok, van der Hulst & Bothun 
1995; van den Hoek et al. 2000), with small-amplitude 
bursts on top of star formation rates that are 
nearly constant when integrated over Gyr timescales 
(Gerritsen & de Blok 1999; Bell et al. 2000; Boissier et al. 
2008; Schombert & McGaugh 2014). The reported lack of 
chemical abundance gradients has been regarded as sup¬ 
porting evidence for star formation histories characterized 
by sporadic episodes, and possibly linked to a mode of 
galactic evolution that does not proceed from the inside out 
(de Blok & van der Hulst 1998; Vorobyov et al. 2009). 

The present work stems from one main motivation. The 
study by de Blok & van der Hulst (1998), which remains to 
date the only investigation of the radial abundance trends in 
LSBGs, reached its conclusion about the lack of detectable 
gradients on the basis of only three galaxies, with a relatively 
small number (4 ^ n ^ 7) of H ll region spectra obtained per 
galaxy. Perhaps more importantly, all three systems have 
an irregular morphological type. Abundance gradients in 
such galaxies are generally found to be very weak or absent 
(e.g. van Zee & Haynes 2006; Croxall et al. 2009), although 
Pilyugin, Grebel & Zinchenko (2015) recently measured sig¬ 
nificant gradients in a number of irregular galaxies. 

Therefore, in order to verily claims, based on the chem¬ 
ical abundance properties, that LSBGs evolve differently 
(aside from the slower evolution) compared to galaxies of 
high surface brightness, it is essential to characterize the 
abundance gradients of LSBGs displaying a spiral morphol¬ 
ogy. Gurrently, multi-object techniques on 8m-class tele¬ 
scopes make the spectroscopy of relatively large samples of 
Hll regions, necessary to obtain reliable estimates of the 
abundance gradients, feasible even with moderate amounts 
of observing time. 

In this paper we present new oxygen abundances of H ll 
regions in 10 late-type LSBGs, with the main goal of deriv¬ 
ing their radial distributions. Our observational material is 
presented in Sect. 2, and we derive nebular chemical abun¬ 
dances using different diagnostics in Sect. 3. We discuss the 
abundance gradients in Sect. 4 and summarize our results 
in Sect. 5. 

2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

LSBGs are often characterized by sparse distributions of 
faint Hll regions, which make it difficult to derive reliable 
chemical abundance gradients. For our study we selected 
galaxies that were shown by previous work to contain a 

^ We adopt 12-|-log(0/H)Q = 8.69 from Asplund et al. (2009). 


sufficient number of Hll regions to alleviate such a diffi¬ 
culty. Our targets were extracted from a variety of sources 
providing imaging and surface photometry information: 
Beijersbergen, de Blok & van der Hulst (1999), Auld et al. 
(2006), de Blok et al. (1995), McGaugh & Bothun (1994) 
and McGaugh, Schombert & Bothun (1995). The galaxies 
selected have R-band central disc surface brightness /rs (0) > 
23 magarcsec”^, a value that is commonly adopted to de¬ 
fine LSBGs (e.g. Impey & Bothun 1997), being more than 
one magnitude fainter than the canonical Freeman (1970) 
value /rs(0) = 21.65 ± 0.3 magarcsec”^ (this is somewhat 
arbitrary, and some authors use, for instance, Ats(O) = 22 
magarcsec“^ as the dividing line between low and high 
surface brightness galaxies, HSBGs; e.g. Burkholder et al. 
2001; Boissier et al. 2003; Liang et al. 2010). We also in¬ 
cluded three targets from the list of extreme late-type spi¬ 
ral galaxies by Matthews & Gallagher (1997). Although the 
value of their central disc surface brightness is not reported 
in their work, these systems have structural and global sur¬ 
face brightness characteristics that are comparable to those 
of the remaining targets in our sample. Their ij,b{0) values 
were estimated from the expression /rs(0) = 25 — r 25 /rd, 
using the values for r 25 and from Lauberts & Valentijn 
(1989), and range between 22.4 and 23.9 magarcsec”^. 

Our final sample of 10 LSBGs comprises systems both 
with and without bulge, and is presented in Table 1, where 
most of the information has been extracted from the Hy- 
perLeda database^ (Makarov et al. 2014). The exceptions 
are reported in the last column of the table. The celes¬ 
tial coordinates of the galactic centers were measured from 
our astrometrically calibrated images (see below). We follow 
the nomenclature from the ESO/Uppsala survey (Lauberts 
1982) and the Uppsala General Catalogue of galaxies (UGC, 
Nilson 1973). We re-classified ESO249-036 (Horologium 
dwarf, part of a pair with ESO 249-035) as an Sm type 
(instead of Im; e.g. Lauberts 1982), based on the patchy, 
arm-like structure shown in the images by Auld et al. (2006, 
Longmore et al. 1982 classified this galaxy as an Sdm). For 
all the remaining galaxies we retain the late-type spiral (Sc 
to Sm) classification given in the sources specified in Table 1. 

We obtained spectra of individual Hll regions in the 
target galaxies with two different instruments and multi¬ 
slit setups: the Very Large Telescope (VLT) equipped with 
the FOcal Reducer and low dispersion Spectrograph (FORS) 
and the Gemini North telescope with the Multi-Object Spec¬ 
trograph (CMOS). The FORS data (collected for the seven 
ESO galaxies) were obtained with the 300V grism using 0'.'9 
slits, while the CMOS data (for the three UGC galaxies) 
were obtained with the B600 grating using l'.'4 slits. 

The slit masks were designed based on narrow-band Ha 
imaging (on- and off-band) secured in advance of the spec¬ 
troscopic runs. For six of the galaxies we used two separate 
masks, in order to increase the number of target H ll regions. 
The number of H ll regions observed per galaxy ranges be¬ 
tween 7 and 21. The total exposure times, reported in Ta¬ 
ble 2, range between 3600 s and 6000 s, and were subdivided 
into a number of sub-exposures varying from three to six. 
The spectral coverage extends from approximately 3650 A 
to 7000 A, with a FWHM spectral resolution of 8 A (FORS) 
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Table 1. Galaxy sample. 


ID 

R.A. (J2000) 

Dec. (J2000) 

Type 

D 

Mb 

PA 

i 

r25 

I’d 

Mfl(O) 

Ref. 


h m s 

o ! ff 


Mpc 


deg 


arcsec 

(kpc) 

mag arcsec“^ 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

ESO 059-009 

07 36 13.08 

-70 42 40.2 

Sd 

16.7 

-17.7 

169 

27 

49 (3.9) 

27 (2.1) 

23.1 

a 

ESO 187-051 

21 07 32.78 

-54 57 09.1 

Sm 

17.8 

-16.8 

11 

52 

40 (3.5) 

26 (2.2) 

23.1 

l,b 

ESO 249-036 

03 59 15.72 

-45 52 20.7 

Sm 

9.7 

-14.9 

156 

15 

35 (1.7) 

37 (1.7) 

24.0 

l,b 

ESO 440-049 

12 05 34.07 

-31 25 25.4 

Sc 

31.2 

-19.0 

62 

42 

45 (6.9) 

18 (2.7) 

23.1 

2,a 

ESO 446-053 

14 21 17.08 

-29 15 47.6 

Sdm 

17.1 

-17.9 

165 

67 

53 (4.4) 

28 (2.3) 

22.9* 

3,c 

ESO 504-025 

11 53 50.15 

-27 20 59.2 

Sd 

24.1 

-17.5 

20 

40 

21 (2.5) 

20 (2.3) 

23.9* 

3,c 

ESO 510-026 

13 57 35.77 

-25 47 27.4 

Scd 

33.4 

-17.8 

130 

37 

31 (5.1) 

13 (2.1) 

22.4* 

3,c 

UGC 1230 

01 45 32.50 

-1-25 31 14.8 

Sc 

57.0 

-17.6 

112 

22 

24 (6.5) 

16 (4.4) 

23.3 

4,5,7,d 

UGC6151 

11 05 56.27 

-1-19 49 34.9 

Sc-Sm 

21.3 

-16.9 

175 

18 

39 (4.1) 

24 (2.5) 

23.2 

6,5,7,d 

UGC 9024 

14 06 40.55 

-1-22 04 12.4 

Sc-Sm 

36.5 

-16.4 

148 

15 

18 (3.3) 

35 (6.2) 

24.5 

5,7,d 


(1) PA, i from Auld et al. (2006). (2) ^s(O) from Beijersbergen et al. (1999), using /rij(O) and the average /^^(O) — IJ-Ft{0) in their 
Table 2. (3) PA, i. Type from Matthews & Gallagher (1997). (4) PA, i from van der Hulst et al. (1993). (5) Type from McGaugh et al. 
(1995). (6) PA, i from Patterson & Thuan (1996). (7) from McGaugh & Bothun (1994). 

Sources for the B-band exponential disc scalelength (a) Beijersbergen et al. (1999) (R-band for ESO440-049); (b) Bell et al. (2000); 
(c) Lauberts &; Valentijn (1989); (d) McGaugh &; Bothun (1994). 

* Not a published value, but derived from the expression /ro(R) = 25 — r^hli'd- 


and 4.8 A (CMOS). Thus all the main nebular emission-line 
diagnostics, from [On] A3727 to [Sii] AA6717,31, are accessi¬ 
ble, and [N ll] A6583 is easily deblended from Ha. 


2.2 Data reduction 

Bias subtraction, flat-field correction and wavelength cal¬ 
ibration were accomplished with the European Southern 
Observatory’s EsoRex pipeline for the FORS data, and 
IRAF® routines in the gemini/gmos package for the CMOS 
data. Flux calibration using standard star spectrophotom¬ 
etry, coaddition of sub-exposures and extraction to one¬ 
dimensional spectra yielded the final data products to be an¬ 
alyzed. Emission line fluxes were measured non-interactively 
with the f itprof s tool in iraf. 

The line fluxes were corrected for interstellar reddening 
by comparing the Ha/H/3 and Hy/H/l ratios to the corre¬ 
sponding case B values at 10"* K, requiring that the two 
line ratios yield a consistent extinction coefficient c(H/l), 
by adjusting the equivalent width of the underlying absorp¬ 
tion component. The reddening-corrected line fluxes for the 
[Oii]A3727, [Oiii]A5007, [Nii]A6583 and [S ii] AA6717,6731 
emission lines, that are used as metallicity diagnostics, are 
reported, in units of /(H/l) = 100, in Table 3 (the printed 
version reports the information for ESO 059-009 only, while 
the complete table is available in the online version of the 
journal). We include the detections of the [Olll] A4363 line, 
which can be used, in combination with [O ill] A5007, to ob¬ 
tain direct estimates of the gas electron temperature. The 
errors quoted include the statistical and flux calibration er¬ 
rors, and the uncertainty in the reddening correction. We 
confirmed that the relative strengths of the emission lines 
are consistent with those expected for photoionized nebu- 


® IRAF is distributed by the National Optical Astronomy Obser¬ 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 


Table 2. Log of the observations. 


Galaxy 

Date 

Exp. time (s) 

ESO 059-009 

2012 Apr 17,18 

6000 (x2 fields) 

ESO 187-051 

2011 June 25,26 

5094 

ESO 249-036 

2011 Mar 04,05 

5094 

ESO 440-049 

2012 Apr 17 

6000 (x2 fields) 

ESO 446-053 

2012 Apr 17,18 

3600 (x2 fields) 

ESO 504-025 

2012 Apr 18 

6000 (x2 fields) 

ESO 510-026 

2012 Apr 17 

6000 

UGC 1230 

2010 Sept 11,12 

3600 (x2 fields) 

UGC 6151 

2011 May 04 

4500 


2012 Mar 20 

4500 

UGC 9024 

2011 June 20,30 

4500 


2012 Feb 26, Apr 15 

4500 


lae, as judged from the appropriate diagnostic diagrams (e.g. 
Baldwin, Phillips & Terlevich 1981). 


3 CHEMICAL ABUNDANCES AND RADIAL 
GRADIENTS 

Measuring nebular chemical abundances that are free 
of large systematic uncertainties remains an unsolved 
problem in astrophysics. The nature of this issue and 
the various attempts to circumvent the difficulties have 
been extensively covered by a large number of au¬ 
thors and will not be repeated here (see, for exam¬ 
ple, Bresolin, Garnett & Kennicutt 2004; Kewley & Ellison 
2008; Lopez-Sanchez et al. 2012). In essence, several emis¬ 
sion line diagnostics and different calibrations of these diag¬ 
nostics have been proposed in the literature, with systematic 
variations on the derived oxygen abundances that reach up 
to 0.7 dex. One additional complication that is especially 
relevant for the current study is the non-monotonic nature 
of some of the diagnostics, which can introduce large un¬ 
certainties when attempting to measure radial abundance 
gradients. 
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Table 3. Reddening-corrected fluxes: ESO 059-009 


ID 

(1) 

R.A. 

(J2000.0) 

(2) 

Dec. 

(J2000.0) 

(3) 

r/r25 

(4) 

[On] 

3727 

(5) 

[Oiii] 

4363 

(6) 

[Oiii] 

5007 

(7) 

[Nil] 

6583 

(8) 

[Sii] 

6731-1-6717 

(9) 

c(H/3) 

(10) 

F(H/3) 

(erg cm“^) 

(11) 

1 

07 

36 

12.4 

-70 

42 

03.50 

0.76 

179 ± 9 


33 ± 2 

71 

± 

4 

76 

± 

3 

0.21 

1.5 

X 

10-16 

2 

07 

36 

08.1 

-70 

42 

09.10 

0.84 

190 ± 10 


37 ± 2 

67 

± 

4 

70 

± 

3 

0.26 

1.1 

X 

10-1® 

3 

07 

36 

10.6 

-70 

42 

10.62 

0.66 

172 ± 9 


45 ± 2 

69 

± 

4 

51 

± 

2 

0.46 

1.5 

X 

10-16 

4 

07 

36 

13.1 

-70 

42 

24.78 

0.32 

116 ± 7 


0 ± 1 

70 

± 

4 

66 

± 

3 

0.43 

7.9 

X 

10-17 

5 

07 

36 

05.9 

-70 

42 

26.15 

0.86 

164 ± 8 


54 ± 2 

66 

± 

4 

59 

± 

2 

0.30 

1.3 

X 

10-1® 

6 

07 

36 

15.5 

-70 

42 

26.62 

0.40 

168 ± 10 


12 ± 1 

72 

± 

4 

52 

± 

2 

0.51 

9.4 

X 

10-17 

7 

07 

36 

15.0 

-70 

42 

38.19 

0.23 

119 ± 6 


18 ± 1 

82 

± 

4 

44 

± 

2 

0.26 

1.6 

X 

10-16 

8 

07 

36 

05.3 

-70 

42 

43.60 

0.90 

194 ± 10 


51 ± 2 

61 

± 

3 

53 

± 

2 

0.45 

2.4 

X 

10-16 

9 

07 

36 

10.5 

-70 

42 

44.25 

0.31 

166 ± 9 


29 ± 2 

96 

± 

5 

83 

± 

3 

0.55 

7.4 

X 

10-17 

10 

07 

36 

14.4 

-70 

42 

46.88 

0.20 

111 ± 6 


15 ± 1 

74 

± 

4 

59 

± 

2 

0.28 

6.4 

X 

10-17 

11 

07 

36 

10.8 

-70 

42 

54.22 

0.40 

188 ± 10 


89 ± 4 

81 

± 

4 

63 

± 

3 

0.35 

5.8 

X 

10-17 

12 

07 

36 

13.1 

-70 

42 

57.41 

0.35 

187 ± 10 


63 ± 3 

71 

± 

4 

40 

± 

2 

0.41 

8.8 

X 

10-17 

13 

07 

36 

16.3 

-70 

43 

07.06 

0.65 

165 ± 8 


14 ± 1 

79 

± 

4 

85 

± 

3 

0.31 

8.7 

X 

10-17 

14 

07 

36 

11.9 

-70 

43 

09.27 

0.62 

96 ± 8 


0 ± 1 

79 

± 

5 

75 

± 

4 

0.12 

5.0 

X 

10-17 

15 

07 

36 

12.5 

-70 

43 

13.57 

0.70 

193 ± 10 


26 ± 1 

84 

± 

5 

84 

± 

3 

0.45 

1.3 

X 

10-16 

16 

07 

36 

16.7 

-70 

43 

36.63 

1.22 

286 ± 18 


70 ± 5 

76 

± 

5 

57 

± 

3 

0.50 

2.6 

X 

10-17 


The complete table, which includes information for the full sample, is available in plain text format in the online version. 
F(H/3) in column (11) is the measured flux of the H/3 emission line. 


A frequently used strategy, and the one we adopt here, 
is to consider a number of different abundance determina¬ 
tion methods, with the intention of identifying the quan¬ 
tities that are largely independent of the diagnostics used. 
For example, radial abundance trends are generally found 
to be qualitatively invariant relative to the choice of abun¬ 
dance diagnostics, but one needs to be aware of the fact 
that different methods can yield different gradient slopes 
(e.g. Bresolin et al. 2009b). In order to compare our results 
with recent studies concerning the statistical properties of 
the gradients of large samples of spiral galaxies, discussed 
below, we have considered the following nebular metallic- 
ity diagnostic methods (as is customary for nebular studies, 
we use the terms metallicity and oxygen abundance inter¬ 
changeably) : 

(i) 03N2 = log[([Oiii]A5007/H,5)/([Nii] A6583/Ha)], with 
the empirical calibration given by Pettini & Pagel (2004), 
which allows a comparison with the study of a large sample 
of spiral galaxies by Sanchez et al. (2014). 

(ii) N202 = [N ii] A6583/[0 ii] A3727, adopting both 

the calibration based on photoionization models by 
Kewley & Dopita (2002, =KD02) and the empirical one 
by Bresolin (2007, =B07). The recent study of abundance 
gradients in local star-forming galaxies by Ho et al. (2015) 
is based on this diagnostic. We have already shown in 
Bresolin et al. (2009a) that the calibrations by KD02 and 
B07 yield abundance gradients having virtually the same 
slopes, despite a large systematic offset. 

(hi) N2 = [N ii] A6583/Hq:, with the calibration by 
Pettini & Pagel (2004). 

(iv) i ?23 = ([Oii]A 3727 -f [Oiii] AA4959,5007)/H,5, adopt¬ 
ing the calibration by McGaugh (1991), in the analytical 
form given by Kobulnicky, Kennicutt & Pizagno (1999). We 
used this diagnostic in order to check whether the metallicity 


gradients derived from the previous diagnostics, all involv¬ 
ing the nitrogen [N ll] A6583 line, are confirmed by consid¬ 
ering only oxygen lines instead. Unfortunately, the use of 
this indicator for abundance gradient studies is complicated 
by the well-known non-monotonic behavior of i? 23 . In order 
to attempt and break the degeneracy we followed other au¬ 
thors in using the [N ll] A6583/[0 ll] A3727 line ratio: objects 
for which [N ll] A6583/[0 ll] A3727 > —1.2 should belong to 
the 7?23 upper branch (e.g. Kewley & Ellison 2008). For our 
sample we found that, in order to obtain monotonic radial 
trends in abundance for most galaxies we need to arbitrarily 
vary this boundary between —1.3 and —1.0. Even so, some 
of the radial trends obtained from 7723 display quite a large 
scatter or even sudden discontinuities, which are unphysi¬ 
cal and largely due to the uncertainty in placing objects on 
the correct branch, and the fact that many of the Hll re¬ 
gions lie in the ‘turnaround’ region of the diagnostic. For our 
purposes this is of secondary importance: we simply wish to 
demonstrate that the radial abundance trends we observe do 
not depend on the selection of diagnostics based on nitrogen 
lines instead of oxygen lines. 

The radial O/H abundance gradients we derive for our 
sample are shown in Fig. 1. The linear least square fits to the 
log(0/H) data obtained from the 03N2 and the N2 meth¬ 
ods yield the same slopes, within the uncertainties, and we 
therefore plot only the results obtained using the former in¬ 
dicator, using green squares for the individual Hll regions, 
and the green line representing the least square fit to these 
data points. For the N202 diagnostic we show the data and 
the corresponding linear hts using both the B07 (red circles 
and line) and the KD02 (yellow triangles and line) calibra¬ 
tions, to illustrate the fact that they yield the same gradient 
slopes. We used the KD02 polynomial expression for the 
N202 index vs. O/H relation calculated for a ionization pa¬ 
rameter q = 2 X lO^cms”^. The N202 diagnostic becomes 
insensitive to oxygen abundance for low O/H values. For 
the KD02 diagnostic we have included in Fig. 1 only data 
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points for which log([N ll] A6583/[0 ll] A3727)> —1.2, corre¬ 
sponding approximately to 12 + log(0/H) > 8.4 in the KD02 
abundance scale. Below such limit the KD02-based O/H val¬ 
ues display large and sudden deviations from the smooth 
radial trends seen in Fig. 1, symptomatic of the breakdown 
of the calibration at low metallicity. We instead extended 
the use of the empirical B07 calibration to a lower limit, 
log([N ll] A6583/[0 ll] A3727)> —1.4, corresponding approxi¬ 
mately to 12 + log(0/H) > 7.8 in the B07 abundance scale. 
Therefore, the number of data points based on the B07 
calibration exceeds that for the KD02 calibration for some 
galaxies (especially ESO 249-036 and UGC 9024, for which 
we are not able to derive a gradient based on the KD02 cal¬ 
ibration; we required a minimum of five data points for a 
fit to the abundance gradients), but we stress again that the 
two calibrations yield slopes that are in very good agreement 
with each other. 

In Fig. 1 we also include, using star symbols, the 
O/H values obtained from the direct method, i.e. based 
on the measurement of the electron temperature using the 
[O III] A4363/A5007 line ratio (this procedure was carried 
out with the nebular package in iraf; see Bresolin et al. 
2009b for details). The weak [Olll]A4363 auroral line was 
measured for 12 Hll regions (six in ESO 446-053). It can 
be seen that the direct abundances are in rough agreement 
with the 03N2 and N202 (B07) diagnostics, as expected, 
since both these calibrations are based on samples of H ll 
regions with direct abundance determinations. Finally, we 
include in Fig. 1 the O/H abundances obtained from R 23 
(cross symbols, without showing the corresponding linear 
fits), but only for those galaxies for which the selection be¬ 
tween the upper and lower branches was unequivocal, based 
on the [N ll] A6583/[Oll] A3727 ratio. In the case of the five 
galaxies for which, over the full radial range covered by the 
H II regions, this line ratio straddles the approximate bound¬ 
ary between the two branches [log([N 11 ] A6583/[0 ll] A3727) 
~ —1.2; e.g. Kewley & Ellison 2008)] we do not show the 
i? 23 -based abundances. 

Table 4 and 5 summarize the slopes and zero points, 
respectively, of the metallicity gradients obtained when the 
nebular O/H values are based on the 03N2 and the N202 
(B07) diagnostics. The slopes are reported, together with 
their errors, in terms of the physical size of the galaxies 
in kpc, of the isophotal radius r 25 and of the B-band ex¬ 
ponential disc scale length rd- We enclose these values in 
brackets when we wish to identify galaxies for which the er¬ 
rors exceed the absolute values of the slope, i.e. when the 
abundance gradient is formally consistent with a flat radial 
distribution. At the bottom of Table 4 we include the av¬ 
erage slopes, with corresponding errors, and their standard 
deviations, calculated including all ten galaxies. 

Before proceeding with the discussion, we make a few 
preliminary remarks, based on the results presented in Fig. 1 
and Table 4: 

(i) a large fraction (50 per cent using the 03N2 diagnostic, 
90 per cent using N202) of the sample galaxies present a 
significant radial abundance gradient The higher number of 
galaxies with gradients that are compatible with being flat 
in the case of 03N2 (5 out of 10) is related to the next point. 
Only one galaxy, ESO 187-051, has a flat gradient according 
to the N202 diagnostic. 


(ii) In Fig. 1 the scatter of the data points corresponding to 
03N2 (the mean rms of the linear fit is 0.095 dex) is approx¬ 
imately two times larger than for N202 (mean rms=0.044). 

(iii) The slopes derived using 03N2 are either comparable 
to (within the errors) or shallower (ESO 249-036, ESO 440- 
049, ESO 510-026) than those obtained from N202. We in¬ 
terpret this difference and (ii) above with the higher depen¬ 
dence of 03N2 on the ionization parameter (e.g. Ho et al. 
2015), and the fact that this diagnostic is constructed with 
both low- and high-excitation emission lines. While it is 
not immediately obvious why 03N2 should yield shallower 
gradients in some galaxies, we also point out that in ab¬ 
solute terms the difference in slope that can be obtained 
between the two diagnostics, while statistically significant, 
is small in our sample of galaxies in relation to the range 
of values found in much larger sample of spirals. Consider¬ 
ing, for example, gradients in units of disc scale lengths the 
average slope from 03N2 is —0.045 ± 0.013 dexr/"^, com¬ 
pared to —0.085 ± 0.010 dexr]/^ from N202 (i.e. the dif¬ 
ference is Ao 3 JV 2 -iV 202 = 0.04 dexr/"^). The slopes found 
by Sanchez et al. (2014) in a sample of 193 spiral galax¬ 
ies (using 03N2) cover a much larger range, from —0.19 to 
+0.04 dexr)/^. Considering isophotal radii instead, the av¬ 
erage slope from 03N2 is —0.065 + 0.018 dexr//^, compared 
to —0.133 + 0.024 dexr^^ from N202 (Ao3jv2 -N202 = 0.07 
dexr^^). For the bulk of the 49 galaxies studied by Ho et al. 
(2015) using N202 the slope ranges from 0 to —0.8 dexr//^. 
Therefore, the small systematic difference resulting from 
the selection of the abundance diagnostic does not hamper 
the comparison with the distribution of abundance gradient 
properties of HSBGs carried out in Sect.4. 

(iv) In those cases where it has been possible to derive a 
reliable radial abundance trend also from R 23 , we find results 
that are compatible either with N202 (all galaxies except 
ESO 249-036) or 03N2 (all galaxies except ESO 440-049). 

For the reasons outlined in points (ii) and (iii) above, 
it is tempting to assign more significance to the results ob¬ 
tained from N202. For completeness, in the following dis¬ 
cussion we will continue to consider gradients obtained from 
both 03N2 and N202. 


4 DISCUSSION 

The first finding of this work has already been introduced in 
the previous section: late-type (spiral) LSBGs do have mea¬ 
surable radial metallicity gradients. We thus do not conhrm 
the earlier result by de Blok & van der Hulst (1998), who 
found no evidence for the presence of such gradients in LS¬ 
BGs. The origin of this discrepancy can be attributed to a 
number of factors that are likely to affect their study: small 
sample size (three galaxies), galaxy morphology (irregulars), 
small number of H ll regions with limited radial coverage, 
and possibly the use of the R 23 diagnostic in the uncertain 
turnaround region. We believe that our result is robust, be¬ 
ing based on a relatively large number of Hll regions per 
galaxy (up to 21), and on deeper spectra than was possible 
in the 1990’s. 

Our study hnally answers the question raised by 
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Figure 1. Radial O/H abundance gradients drawn in terms of the isophotal radius r 25 (lower scale of each panel) and of the physical 
radius in kpc (upper scale) from different diagnostics (see symbol explanation in the top left panel). In addition, star symbols are used for 
abundances obtained with the direct method (when available), and crosses for abundances from R 23 (only shown when branch selection 
is unambiguous). The lines represent least square fits to the abundance data. See text for explanations. 
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Table 4. Abundance gradient slopes. 


ID 

dex kpc ^ 

cr 

03N2 
dex r 2 ^ a 

dex ^ 

cr 

dex kpc ^ 

cr 

N202 

dex r 2 ^ a 

dex ^ 

cr 

ESO 059-009 

-0.042 

0.020 

-0.165 

0.080 

-0.090 

0.044 

-0.040 

0.010 

-0.156 

0.039 

-0.085 

0.021 

ESO 187-051 

(-1-0.005 

0.023 

-1-0.019 

0.079 

+0.012 

0.051) 

(-0.004 

0.010 

-0.014 

0.035 

-0.009 

0.022) 

ESO 249-036 

(-1-0.001 

0.040 

-1-0.001 

0.066 

+0.001 

0.070) 

-0.065 

0.018 

-0.108 

0.030 

-0.114 

0.031 

ESO 440-049 

-0.013 

0.011 

-0.092 

0.076 

-0.036 

0.030 

-0.037 

0.005 

-0.256 

0.035 

-0.101 

0.014 

ESO 446-053 

(-0.006 

0.031 

-0.028 

0.138 

-0.015 

0.073) 

-0.023 

0.017 

-0.102 

0.074 

-0.054 

0.039 

ESO 504-025 

(-0.025 

0.028 

-0.063 

0.069 

-0.060 

0.065) 

-0.034 

0.011 

-0.086 

0.028 

-0.081 

0.026 

ESO 510-026 

(-0.009 

0.017 

-0.047 

0.087 

-0.019 

0.035) 

-0.047 

0.009 

-0.239 

0.047 

-0.097 

0.019 

UGC 1230 

-0.022 

0.007 

-0.145 

0.043 

-0.099 

0.029 

-0.023 

0.004 

-0.149 

0.026 

-0.101 

0.018 

UGC 6151 

-0.020 

0.017 

-0.082 

0.069 

-0.050 

0.042 

-0.039 

0.006 

-0.161 

0.026 

-0.098 

0.016 

UGC 9024 

-0.015 

0.006 

-0.049 

0.020 

-0.092 

0.038 

-0.018 

0.004 

-0.059 

0.012 

-0.112 

0.022 

Average 

-0.015 

0.004 

-0.065 

0.018 

-0.045 

0.013 

-0.033 

0.005 

-0.133 

0.024 

-0.085 

0.010 

Standard dev. 

0.013 


0.055 


0.038 


0.016 


0.072 


0.030 



The brackets are used to identify regressions for which the slope is flat within the quoted errors. 


Table 5. Abundance gradient zero points. 


ID 

03N2 

12+log(0/H) 

N202 

12+log(0/H) (T 

ESO 059-009 

8.796 

0.054 

8.607 

0.025 

ESO 187-051 

8.256 

0.076 

8.052 

0.033 

ESO 249-036 

8.181 

0.071 

8.046 

0.032 

ESO 440-049 

8.543 

0.062 

8.434 

0.029 

ESO 446-053 

8.227 

0.075 

8.122 

0.040 

ESO 504-025 

8.439 

0.105 

8.263 

0.043 

ESO 510-026 

8.355 

0.080 

8.255 

0.043 

UGC 1230 

8.428 

0.061 

8.206 

0.035 

UGC 6151 

8.339 

0.052 

8.098 

0.019 

UGC 9024 

8.308 

0.055 

8.049 

0.029 


Edmunds (1999), who speculated that the chemical evolu¬ 
tion of low and high surface brightness galaxies should pro¬ 
ceed in a similar fashion, and wondered whether LSBGs and 
HSBGs ‘might not show similar radial abundance gradients, 
since organization of star formation by spiral structure is 
a good candidate mechanism for gradient generation’. In¬ 
deed our observations show that, as in the case of HSBG 
spirals, metallicity gradients are a common characteristic 
of spiral LSBGs, implying that models of inside-out disc 
formation can explain the chemical abundance properties 
of these galaxies, too. It is not necessary to invoke dras¬ 
tically different modes of disc buildup for LSBGs, as sug¬ 
gested in the literature in order to account for the apparent 
lack of abundance gradients (de Blok & van der Hulst 1998; 
Vorobyov et al. 2009). 

The question we wish to address next is whether 
the gradients observed in LSBGs are quantitatively com¬ 
parable to those of HSBGs. For this purpose we com¬ 
pare the slopes obtained in Section 3 with those mea¬ 
sured from large samples of HSBGs recently published 
by Sanchez et al. (2014) and Ho et al. (2015). Both 
studies confirm, with higher statistical signihcance, ear¬ 
lier Endings obtained from smaller samples of galax¬ 
ies and Hll regions (e.g. Vila-Gostas & Edmunds 1992; 
Zaritsky, Kennicutt & Huchra 1994; Garnett et al. 1997), 
namely that spirals are characterized by similar abundance 


gradients, irrespective of morphological type and other 
global properties, when normalized to either the exponen¬ 
tial disc scale length or the isophotal radius. 

a) Sanchez et al. (2014) analyzed the abundance gradient 
properties of 306 galaxies observed as part of the Calar Alto 
Legacy Integral Field Area survey (CALIFA, Sanchez et al. 
2012). Ghemical abundances could be derived for at least 
four Hll regions per system in 193 galaxies, allowing a de¬ 
termination of their radial abundance gradients. These au¬ 
thors established the existence of a characteristic slope of 
the radial metallicity (O/H) gradient, that is independent 
of the properties of non-interacting, isolated galaxies, such 
as morphological type, presence of bars, and stellar mass. 
Considering only their subsample of 146 isolated galaxies, 
showing no evidence for ongoing interactions, the character¬ 
istic gradient, expressed in terms of the disc eEective radius 
re, has a value a = —0.11 dexrj^, with a standard devia¬ 
tion of 0.08 dexr“^. The normalized histogram of the slope 
distribution, extracted from Fig. 8 in Sanchez et al. (2014), 
is shown in the top panel of Fig. 2 by the black, dot-dashed 
line. 

We can make a direct comparison with the average value 
we obtained for LSBGs using the 03N2 diagnostic, which 
was also adopted by Sanchez et al. (2014): ar^ = —0.045 ± 
0.038 dexr^^ (mean ± standard deviation), which translates 
into are = —0.076 ± 0.064 dexr“^ (accounting for the fact 
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Figure 2. (Top) Histogram of the slope distribution for the LS- 
BGs (red continuous line) and for the sample by Sanchez et al. 
(2014, black dot-dashed line). Units for the slope are dexr^T^. 
(Bottom) Histogram of the slope distribution for the LSBGs (red 
continuous line) and for the sample by Ho et al. (2015, black dot- 
dashed line). Units for the slope are dexr^^. In both panels mean 
values (diamonds) and standard deviations (horizontal bars) are 
indicated for the different samples. 

that Te = 1.678 rji). The distribution of the slope values we 
obtain for the LSBGs is shown in the top panel of Fig. 2 by 
the red, continuous line. It can be seen that the mean slopes 
for the HSBGs and LSBGs are indistinguishable, and that 
the two distributions are also compatible with each other, as 
confirmed by the t-test and the Kolgomorov-Smirnov test, 
respectively. 

We note that, while the GALIFA sample analyzed by 
Sanchez et al. (2014) is complete only at the bright end, the 
luminosity of the 146 non-interacting galaxies that consti¬ 
tute our comparison sample extends down to Mg ~ —17, 
thus overlapping with the luminosity range of the LSBGs in 
our sample. Moreover, the analysis by Sanchez et al. (2014) 
excluded trends of abundance gradient slopes with galaxy 
luminosity. 

b) Ho et al. (2015) presented the metallicity gradient prop¬ 
erties of 49 spiral galaxies with absolute magnitudes —16 > 
Mb > —22 in the N202 (KD02) abundance scale. We note 
the good overlap in absolute magnitude with our LSBG sam¬ 
ple. In the previous section we showed that when measuring 
gradient slopes, in alternative to the KD02 calibration, we 
can use the B07 calibration, with the results summarized 
in Table 4. Ho et al. (2015) found evidence for a common 
gradient among their galaxies, when the slope is expressed 
in units of the isophotal radius r 25 . For the fainter half of 
their sample {Mb > —20.1), to which we compare our galax¬ 
ies, they found a mean slope aras = —0.34 dex with a 
standard deviation of 0.17 (for the full sample the slope is 
similar: —0.39 dexr^g^, with a standard deviation of 0.18). 
This compares to an average value ar 25 ~ —0.133 dexr^g^ 
for our LSBG sample, with a standard deviation of 0.072. 



Figure 3. Slope of the abundance gradients, in dexkpc“^, as a 
function of the H-band exponential disc scale length for our LS¬ 
BGs (red full dots) and a sample of HSBGs (open squares; slopes 
from Pilyugin et al. 2014a and disc lengths from Pilyugin et al. 
2014b). 

The bottom panel of Fig. 2 shows the distribution of the 
gradient slopes for the Ho et al. (2015) sample (black, dot- 
dashed line) and the LSBGs (red, continuous line). 

The t-test indicates that in this case the difference be¬ 
tween the average slopes of HSBGs and LSBGs is statisti¬ 
cally significant {p = 0.0004). A similar comparison with the 
HSBG sample presented by Pilyugin et al. (2014a), exclud¬ 
ing galaxies in pairs or mergers, as discussed by Ho et al. 
2015, yields the same conclusion. 

4.1 Interpretation of the observed trends: clues 

on the common abundance gradients of spiral 
galaxies 

We show in Fig. 3 the slope of the abundance gradients, now 
expressed in dexkpc“^, as a function of the reciprocal of the 
B-band exponential disc scale length 1/rd for our LSBGs 
(red full dots) and a sample of HSBGs for which the relevant 
information is readily available (open squares; slopes from 
Pilyugin et al. 2014a and disc lengths from Pilyugin et al. 
2014b). The plot illustrates how LSBGs follow the same 
trend observed for HSBGs, namely that the gradient slope in 
dexkpc”^ anticorrelates with galaxy size: large discs display 
shallower gradients compared to small discs. As suggested by 
Prantzos & Boissier (2000), who discussed a similar plot in 
their modeling of spiral galaxies as a function of the rota¬ 
tional velocity and spin parameter of their haloes, this trend 
describes a fundamental property of galaxy discs, since 1 /rd 
measures the radial decrease of the surface brightness with 
radius in mag arcsec”^ kpc“^ [this follows from the equa¬ 
tion for the exponential decrease of the surface brightness 
E(r) = E(0) exp(—r/rd) and the corresponding equation 
when using magnitudes /r(r) = po + 1.086r/r^]. Our plot 
in Fig. 3 suggests that the underlying cause for the rate 
of radial decrease of the metallicity per kpc in discs is the 
same for all galaxies, regardless of their surface brightness. 
We argue that this is simply the rate of variation of the 
mass surface density with radius (also described by 1/rd), 
through the local mass surface density-metallicity rela¬ 
tion Em — Z (Edmunds & Pagel 1984; Rosales-Ortega et al. 
2012; Moran et al. 2012), or equivalently the local stellar 
surface brightness-metallicity relation El — ^ (Ryder 1995; 
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Pilyugin et al. 2014b). Smaller galaxies have a steeper slope 
(in dexkpc”^) in the radial decrease of both the surface mass 
density (or surface luminosity) and metallicity. A linear least 
square fit to the data in Fig. 3 yields the following: 

dlog{0/H) ^ _Q ^Q2 (±0.014)— ± 0.002 (±0.005) (1) 

dr Td 

with distances expressed in kpc. Recalling that we can write 
l/vd = A/i/(1.086Ar) using the variation of the surface 
brightness between two points separated by a linear distance 
Ar, we can derive from Eq. (1) a simple relation between the 
change oi O/H and the corresponding variation in surface 
luminosity log E: 

Alog(C>///) = 0.235 (±0.032) A log (2) 

From the same HSBG galaxies included in Fig. 3 
Pilyugin et al. (2014b) determined that Alog(0/Ff) = 
0.286 (±0.090) A log E_l^ for the centers of spiral discs, 
which agrees with Eq. (2) within the errors. Relations with 
comparable coefficients, valid across the discs of spirals, were 
found by Ryder (1995), but for different photometric pass- 
bands (e.g. 0.237 ± 0.085 in the /-band). 

If the underlying mechanism is the same why, then, is 
the mean slope of the gradient of LSBGs shallower than 
for HSBGs, when expressed in terms of the isophotal radius 
r 25 , but similar when referred to the exponential disc scale 
length, as shown above? Our interpretation, outlined below, 
is that the former result is due to a sort of selection effect, 
dependent on the disc central surface brightness. Our argu¬ 
ment also explains the existence of a common abundance 
gradient slope for HSBGs, as confirmed by Sanchez et al. 
(2014) and others. 

Let us consider abundance gradients expressed in terms 
of the isophotal radius, r 25 . For a typical HSBG, we adopt 
the central surface brightness as given by the canonical 
Freeman (1970) value of Ms(0) = 21.65 magarcsec”^. 
Across the isophotal radius, the surface brightness varies 
by A/ 4 jg 3 Q = 25 — 21.65 = 3.35 magarcsec”^. For a typ¬ 
ical LSBG, we adopt instead Ub{0) = 23.2 magarcsec”^ 
(the median value in Table 1, but this choice is some¬ 
what arbitrary), and the variation over the isophotal ra¬ 
dius is reduced to A/j^sbc ~ magarcsec”^. Convert¬ 
ing surface brightness to surface luminosity El (e.g. mea¬ 
sured in Lqpc~^) the corresponding radial variations are 
(AlogEL)HSBG = 1-34 and (AlogEL)LSBG = 0-72- We as¬ 
sume now that the linear form of the local relation between 
log(0/H) and log El applies to all galaxies. We avoid the 
use of the local Em — Z relation, which requires an assump¬ 
tion about the mass-to-light ratio, and although these sur¬ 
face density-metallicity relations are better fit with second 
order polynomials (Rosales-Ortega et al. 2012) or using ad¬ 
ditional parameters (e.g. disc scale length or morphology, 
Pilyugin et al. 2014b), the intent of our estimate is simply 
to reproduce the essential trends. 

Using Eq. (2) we can now estimate the variation 
of the O/H ratio across the isophotal radius r 25 to be 
Alog(0///)HSBQ = 0.31 dex and Alog{0/H)^g^^ = 0.17 
dex. These estimates compare well with the actual average 
measurements discussed above (0.34 and 0.13 dex, respec¬ 
tively), and the small differences can easily be attributed 


both to the approximate nature of our calculations (includ¬ 
ing the choice of the central surface brightness for LSBGs) 
and the relatively large scatter in the observational data. 
We can conclude that LSBGs have a shallower mean abun¬ 
dance gradient in comparison to HSBGs, when expressed in 
terms of the isophotal radius, simply as a result of the re¬ 
duced range in surface brightness measured between their 
centers and their isophotal radii. This effect is naturally 
not observed when the gradients are expressed in terms of 
the exponential disc scale length because over such a length 
the variations in surface brightness, surface luminosity and 
surface mass density are independent of the central surface 
brightness (e.g. the surface brightness decreases by a con¬ 
stant 1.086 magnitudes arcsec”^ over this length). 

A far-reaching corollary to these simple considerations 
is that the common abundance gradient slopes observed in 
HSBGs, expressed either in terms of disc scale lengths or 
isophotal radii (as found by Sanchez et al. 2014, Ho et al. 
2015 and others), and their relatively small dispersions, de¬ 
rive naturally from the surface mass density (or luminosity)- 
metallicity relation (see also the discussion in Sanchez et al. 
2014) and, in the case of the gradients per isophotal radii, the 
narrow range in central surface brightness spanned by the 
HSBGs typically selected for chemical abundance analysis. 
With a dispersion of 0.3 magarcsec”^ in ^b(O) (Freeman 
1970) and a dispersion of 0.11 dex in O/H coming from 
the empirical "Zl — Z relation (Pilyugin et al. 2014b) we 
roughly estimate a dispersion in the gradient slope of ~ 0.12 
dexr^g^, which compares well with a — 0.12 dexrjg^ and 
a = 0.18 dexr^g^ reported by Sanchez et al. (2014) and 
Ho et al. (2015), respectively. These figures can be improved 
with appropriate modeling of the uncertainties involved, but 
this lies beyond the scope of our study. Considering scale 
lengths instead, the mean value found in the ^-band by 
Fathi et al. (2010) for a large sample of galaxies from the 
Sloan Digital Sky Survey is 3.85 ± 2.10 kpc. Using Eq. (1) 
we derive a corresponding gradient of a = —0.15 dexr^^, 
with a standard deviation of 0.08 dexr“^. There is a system¬ 
atic offset in comparison with the average slope measured 
by Sanchez et al. (2014), but the standard deviation is the 
same. 

Lastly, we conclude that the common abundance gra¬ 
dient measured in dex expresses a more fundamental 
property of spiral discs than the common gradient in dex 
r/g^, being linked to the underlying physics (the existence 
of a local mass surface density-metallicity relation and its 
relative variations with radius) in a simpler way, without de¬ 
pendence on structural properties such as the actual value 
of the central mass surface density. 


5 SUMMARY 

Our observations show that LSBGs have measurable radial 
chemical abundance gradients, qualitatively similar to those 
present in the more extensively studied high surface bright¬ 
ness galaxies. This result contrasts with the previously held 
notion that LSBGs do not possess such gradients. In ad¬ 
dition, we found that the mean LSBG gradient, when ex¬ 
pressed in terms of the isophotal radius, is significantly shal¬ 
lower than for the average HSBG. Since both LSBGs and 
HSBGs follow the same relation between the radial gradients 
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of both metallicity (in dexkpc“^) and surface brightness (in 
magnitudes arcsec”^ kpc“^), we can explain this result sim¬ 
ply adopting a common surface luminosity-metallicity rela¬ 
tion (reflection of the more fundamental relation involving 
the mass surface density) and considering the range in sur¬ 
face brightness across the optical discs. The common gradi¬ 
ent in dex 1-25^ observed in HSBGs derives from the similar¬ 
ity of their surface brightness (or equivalently mass surface 
density) gradients and central surface brightness. The lat¬ 
ter conclusion does not contradict alternative explanations 
for the existence of this common abundance gradient. For 
example, Ho et al. (2015) based their interpretation on the 
similarity of the radial gas and stellar surface density profiles 
among spiral discs. Our reformulation, based on the struc¬ 
tural similarity and the surface brightness properties, yields 
a concordant, equivalent interpretation. On the other hand, 
when normalizing the gradients to the exponential disc scale 
lengths we directly expose the underlying physical mecha¬ 
nism, i.e. the mass surface density-metallicity relation, and 
the homology with which it varies with galactocentric dis¬ 
tance among galaxies. 
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